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The coupling of modes in two parallel dielectric waveguides is studied. 
The individual waveguides are assumed to be asymmetric and unlike each 
other. If the individual waveguides support modes with nearly equal 
propagation constants /3 2 and /3 4 = /3 2 + IB A, then the double waveguide 
system will support two new modes with propagation constants /3_ = /3 2 — 8 
and /3+ = /3 4 -f- 5. The shift 8 is related to A and to the shift 8 which would 
occur if the original modes were degenerate; 8 is expressed in terms of the 
parameters describing the asymmetric double waveguide system. The field 
distributions of the new modes are approximately even and odd combina- 
tions of those of the original modes in the isolated waveguides; the relative 
amplitudes with which they are combined depend upon the amount of 
mismatching A. As the modes travel down the waveguide system, they 
partially cancel and add, thus transferring power. A power transfer ratio 
F is defined and is shown to decrease rapidly as A/5 increases. The beat 
length L depends upon both 8 and A/8; it also decreases as A/6 increases. 
A numerical example is given to illustrate the effects of mismatching and 
to demonstrate the feasibility of constructing a mode-coupling device. 
Possibilities of tuning the device to reduce mismatching are discussed. 

I. INTRODUCTION 

Coupling of degenerate modes of parallel optical waveguides has 
been discussed by Kapany 1 and, to a greater extent, by Marcuse. 2 
Such coupling is of particular interest in the field of fiber optics, since 
it may cause undesirable crosstalk between adjacent optical fibers used 
for light transmission. Marcuse 2 has applied the theory of degenerate 
mode coupling to the problem of crosstalk between cladded optical 
fibers embedded in a lossy medium and between cladded dielectric 
slab waveguides. The fabrication of devices which would actually 
take advantage of mode coupling, such as for light switching, modula- 
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tion, or power transferral, 3 is fraught with practical difficulties, since 
the specification of physical parameters must necessarily be stringent. 
These difficulties require us to view the theory of optical waveguide 
coupling from a new vantage point. 

Let us first sketch briefly what is known. If two optical waveguides 
each have a mode with the same propagation constant /3, then when 
the two waveguides are placed parallel to each other, the double wave- 
guide system supports two new modes whose propagation constants 
are /3+ = /3 + 5 and |8_ = /3 — 5. These two modes are approximately 
symmetric and anti-symmetric combinations of the original modes in 
the isolated waveguides. The shift in propagation constant, 5, is 
related to the coupling coefficients involved in a description of the 
modes by means of general coupled line equations. It can also be ex- 
pressed via a perturbation treatment of Maxwell's equations. Since 
the superimposed modes travel down the double waveguide system at 
different phase velocities, they alternately add and cancel. If the wave- 
guides are lossless, power is transferred back and forth over a beat 
length L = 71-/ (25). On the other hand, Marcuse shows that, if the 
waveguides are lossy, they tend to equalize the power they carry, 
provided the modes travel far enough. A lossy external medium also 
causes mode loss. Marcuse further states that only degenerate modes 
exchange a significant amount of power if their coupling mechanism 
is independent of length. 

With this abbreviated version of the present theory in mind, we see 
several criteria which a mode coupling device should satisfy: (i) the 
core and cladding of each waveguide should be lossless, (ii) the medium 
external to the waveguides should be lossless, and (Hi) the two wave- 
guides should have a degenerate mode. (There are also other criteria, 
such as that the waveguide walls be free from imperfections, but they 
are not discussed.) 

The first criterion is an important one and certainly merits further 
study. In this paper, however, we avoid the issue by assuming that the 
device we fabricate has lossless waveguides. A subtler way to put this 
is to say that the device is short enough that losses can be ignored. 

The second criterion is satisfied by assuming that the claddings of 
the two waveguides are contiguous and that there is no medium ex- 
ternal to them. Instead of thinking in terms of two optical fibers, we 
consider two dielectric slab waveguides placed next to each other. 
Fabrication would be similar to that currently used in the production 
of double heterostructure lasers and modulators. 4,5 Each waveguide 
will consist of a slab of high refractive index surrounded by two slabs 
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of lower index. Since the double waveguide device will have a central 
slab common to both waveguides, the device can be modelled by a 
5-dielectric-slab model. 

The third criterion, that the two waveguides have a degenerate mode, 
motivates our present study. In practice, it is very difficult to fabricate 
a device with degenerate modes. It is therefore quite important to 
know how well the device will operate if the propagation constants for 
the modes are slightly mismatched. We study the effect of mismatching 
on the beat length and on the capability of the device to transfer power. 
We also discuss methods of tuning the device after it is fabricated. The 
tuning could be used to match the propagation constants more closely. 
It might also be used dynamically, thus offering the possibility of 
utilizing the double waveguide system as a light switch or a modulator. 

II. FORMULATION 

We adopt the standard slab model of an absorptionless dielectric 
medium. The optical dielectric K(x), i.e., the square of the refractive 
index, is assumed to vary only with x and to take the piecewise con- 
stant form shown in Fig. 1. If the waves are assumed to travel in the 
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Fig. 1 — The optical dielectric profile K{x) for the five-slab model. 
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z-direction with propagation constant /3, then the electric and magnetic 
fields are independent of y and can be expressed as 

E = e(x) exp i(ut — @z), 
H = h{x) exp i(u)t — 0z), 

where o> is the angular frequency of the light and t is the time. Both 
TE and TM modes exist. It follows from Maxwell's equations that the 
electric field e y (x) of a TE mode is described by 

|£ + Zk*K(x) - p-]e v = 0, (1) 

and that the magnetic field h y {x) of a TM mode is described by 

It is required that 



h. 



de y 
dx ' 

1 dh y 
K(x)~dx~ 



be continuous. Since K(x) is piecewise constant, solution of eqs. (1) 
and (2) subject to the above conditions is straightforward. The 
solution of (1) is 

e v (x) = A exp p& x < 

= A[(pi/p 2 ) sin p 2 x + cos p 2 x~\ < x < 2w 2 
= AC 2 [1 + (p 1 /p 2 )7 1 2 ][-Zsinhp 3 (x - 2w 2 ) 

+ cosh p 3 (x — 2w 2 )~\ 2w 2 < x < 2(w 2 + w 3 ) 
= AC 2 C 3 tl + (pi/p«)r,][l - XT{] 

X [(p 3 /p 4 )Fsinp 4 (x - 2w 2 - 2w z ) + cosp 4 (x - 2w 2 - 2u> B )] 
2(to s + w 3 ) < re < 2(w 2 + w 3 + w t ) 
= AC 2 C 3 C 4 [1 + (pi/p 2 )7 T 2 ][l - Xr,] 

X [1 + (piMYTt] exp p 6 (2w 2 + 2w 3 + 2w 4 - a) 

2(w 2 + ^3 + w A ) < x, (3) 

with 

PiiP) = {p - k*K& i = 1, 3, 5, 

p,-03) = {WKi - /3 2 )* i = 2, 4, 

C»(0) = cos2p 2 w 2 , T»{0) = tan2p 2 w 2 , 

C«(/3) = cosh2p 3 u; 3 , r 8 Q3) = tanh 2p 8 w 8 , (5) 

CM = cos2p 4 w 4 , T A (P) = tan2p 4 w 4 , 
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(4) 



'» -strife]- 

The amplitude .A is arbitrary. Equation (3) satisfies the continuity 
condition on e y (a;) everywhere, and that on de y /dx at all but the point 
x = 2{w 2 + W3). The continuity condition at this point leads to the 
eigenvalue equation 

t (8) - X W + gg> ™ 

which determines the values of the propagation constant 8 for which 
discrete modes can exist. 

For TM modes, the analogous equations are formed by replacing 
Pi by pi and Wi by w„ where 

KipS) = <P - VKtf i = 1, 3, 5, 
KMB) = (k 2 Ki - /3 2 )l i = 2, 4, 
Wi = KiWi i = 2, 3, 4, 

and /3 denotes the propagation constant for a TM mode. For simplicity 
of exposition, we have mainly confined our analysis to that of TE 
modes. It should be clear how to do the corresponding analysis for 
TM modes. 

We remark that the above analysis is quite general and makes no 
assumptions about the relative heights or widths involved in the 
dielectric profile K(x) sketched in Fig. 1. By making appropriate 
choices of the parameters, we could deduce from eqs. (3) to (8) the 
corresponding equations for a single asymmetric or symmetric wave- 
guide, for example. Two cases which interest us particularly are: (i) 
K 3 = Ki = K 5 with K 2 > K x and K 2 > K 3 , and (it) K x = K 2 = K 3 
with Ki > Kz and Ki > K 5 . Each of these models an isolated wave- 
guide. We call the first of these (with high dielectric region K 2 ) guide 
II and the other (with high dielectric region K t ) guide IV. The eigen- 
value equation (8) then reduces to 

X(B) = 1 guide II, 

Y(8) = 1 guide IV. 

These may appear more familiar to the reader when cast in the standard 

MODE COUPLING IN WAVEGUIDES 721 



form for a single asymmetric guide 6 

4. o (Pi/Pi) + (P3/P2) • , TT , ft v 

tan 2wm = 1 - Wp*Hp*M gmde "' (9) 

tan 2w iPi = WP4)+( P 5/r4) ide Iy (1Q) 

1 - (P3/P4) (P5/P4) 

For the model of two adjacent waveguides, we place guides II 
and IV adjacent to each other as illustrated in Fig. 1, with 
K 2 > max [Ki, K 3 ], K t > max \K 3 , K 5 ), and u>< > (i = 1, 2, 3, 4). 

It is also a relatively straightforward procedure to write down 
precisely how many modes can exist with a given dielectric profile 
K(x). 7 Although we omit such expressions here, we do comment that, 
just as a single asymmetric guide may not be able to support a propa- 
gating mode, so also an asymmetric double waveguide structure is 
not always capable of mode propagation. If Ki = K 3 = K 6) though, 
so that the structure is composed of two parallel symmetric (but not 
necessarily identical) waveguides, then there is always at least one 
mode. 

III. NEARLY DEGENERATE MODES 

Let 02 and 4 denote solutions of X(@) = 1 and F(0) = 1, respec- 
tively ; 2 and /3 4 , then, are propagation constants for modes in guides 
II and IV if the guides were isolated from each other. We need make 
no assumption about the order of each mode. In practice, though, both 
propagation constants are likely to be associated with zeroth order 
modes. For definiteness in notation, we assume that /3 4 ^ 02 and write* 

/3 4 - 02 = 2A. 

We now assume that A is "small," i.e., that the two modes are nearly 
degenerate. This assumption, which is fundamental to the remainder 
of the analysis, is stated more explicitly later [eq. (18)]. Frequently, 
it does not matter (to the order of approximation used) whether 2 or 
04 is used in the evaluation of an expression. In such instances, it 
sometimes is helpful to use the notation 0o(— 02 — 04) • 

Our task is now to determine values of which satisfy (8). Our 
experience with the degenerate case (A = 0) leads us to expect that 



* In the numerical example of Section VI, we relax this notation to read 
1 04 — 0i I = 2A, where it is not known a priori whether 0* or 04 is larger. This 
should not be confusing when taken in context. 
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there will be two solutions 0+ and 0- close to O . A study of the coupled 
line equations 8 for the two modes would demonstrate that 0+ = 04 + 8 
and 0- = 0i — 8, where 8 is expressed in terms of the (unknown) cou- 
pling coefficients. 6 We prefer to attack (8) directly; we shall verify the 
expressions for 0+ and 0-, prove that 8 > 0, and give an explicit formula 
for I 

We first show that, if /3 2 and 0a are close enough together, then (8) 
has no solution such that 02 ^ ^ 0*. Since X(0 2 ) = 1 and 
X'(0)£O for all 0, we know that if (0-0 2 )/02«l, then 0<XQ3) £1. 
Similarly, Y(0) ^ 1. Thus, 

X(0) + Y(0) 1 +X(/3)F-i(/3) > j 



l+X(/3)F(/3) Z(/3) + y-»(/5) 

But ZV/S) = tanh 2w 3 p3 < 1 for all 0, so (8) is not satisfied. 

Next, suppose (8) has a solution + = 0i + 5, with 5 > 0. Then if 
5 is small, we know that < X(/3+) < 1 and < Y(0+) < 1, so (8) 
may be written as 



T(/3) ■ tanh w 3 p 3 = tanh [^tanh -1 X + tanh" 1 F)] 

X[l + (1 - F 2 )*] + ^[1 + (1 - X 2 )*] 
XK + [1 + (1 - X 2 )»][l + (1 - F 2 )*]' 

We have 

X(/3+) = X(0, + 2A + 5) = 1 + (2A + 8)X'(0 O ), 
Y(0 + ) = F(/3 4 +5) = l + 5F'(/3o). 



(ID 



(12) 



If we substitute (12) in (11) and perform a perturbation analysis 
under the two assumptions, 

[5(2A + 8)X'(0 o )Y'(0 o )y « 1, (13) 

[-(2A + 8)X'(0 O )J> + t-~8Y'(0 o )J> « 1, (14) 

we find 

TOo) = 1 - [-(2A + 8)X'(0 o )y\:-8Y'(0o)y, 

so that 

8 = -A+ (A 2 + 5 2 )', (15) 



when; 



_ 1 - T(0 O ) 
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On the other hand, if we suppose that (8) has a solution £_ = /3 2 — £, 
with $ > 0, then X(/3_) > 1, F(0_) > 1, and (8) becomes 

5T(/9) ■ tanh iu 3 p 3 = tanh OKtanh -1 X~ l + tann" 1 F" 1 )] 

X + (X 2 - 1)* + Y + (F 2 - 1)* 
= 1 + [X + (X 2 - 1)»][F + (F 2 - 1)*]" 

By a procedure very much like that used to determine 5 = /3+ — /3 4 , 
we find that 8 = 8, as anticipated. Here, the roles of X' and Y' must 
be interchanged in (14). 

The effect, then, of placing guides II and IV next to each other is to 
shift their (isolated) propagation constants /3 2 and /3 4 symmetrically 
outward by 8 to /3_ and f3+. The physical meaning of 5 in (16) is clear: 
it is the magnitude of the shift which would occur if guides II and IV 
had degenerate modes (A = 0). We shall call 8 the "degenerate shift." 

Let us consider assumptions (13) and (14) in more detail. By means 
of (15) and (16), (13) becomes 

1 - r(j8 ) = 1 - tanh w a p 3 « 1. (17) 

This, then, is essentially a restriction on the separation between the 
two waveguides. If they are too close, our approximations will break 
down. Assumption (14) and its counterpart with the roles X' and Y' 
reversed are, by (15), satisfied if 

[A + (A 2 + «»)»]*{[-X'05.)]» + [-F'OSo)]*} « I- (18) 

This tells us how large A can get without invalidating the approxi- 
mations. 

By using eqs. (4), (6), and (7) cleverly, we see that the expressions 
for X'(/3o) and Y'(@ ) reduce to the simple forms 

X'03o) = ~^ P IJ~ Vl) {2 Pl w 2 + 1 + (pi/p01 d9) 

P1P2P3 

F'(0o) = ~ go(pl g + Vl) O5W4 + 1 + (p./pi)]. (20) 

Thus, by (16), 

p 2 p3P4[l —tanh wzp{] 



8 = 



. B (pi+pJ)(p!+pJ)[2pi«> J + H-g J 



2p 5 i<;4+l+- 

P3j 



(21) 
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If we take the case of two identical symmetric waveguides 
(pi = Pa = Ps, P2 = p*, w 2 = Wt) and use the approximation 

tanh-u> 3 pi — 1 — 2 exp ( — 2w 3 pi), 

then (21) reduces to 

& ^ pfplexp (-2w 8 pi) .„„. 

/3o(p? + p 2 2 )(l+pi^ 2 ) W 

which is in agreement with results of Marcuse. 2 For TM modes, we 
arrive at 

B+ = & + I, B- = fa - I 2 = ft - &, 

where 

2 = -2 + (A 2 + £*)*, 

* _ 1 — tanh w 3 p 3 

( K\p\ + K\pl \ h ( Klpl + Klpl \ 1 
+ M JTRi ) + h\ W+Tl )]' 

, ^ 2 / jqg + ggg \ , fr / gfgj + gag \ i 

+ M e + a J + sv « + a )y 

IV. A LOOK AT THE MODES 

We now discuss what the modes e± (x) associated with /3 ± are like in 
guides II and IV. The expressions for e±{x) are given by (3). Both the 
shapes of the modes and their relative amplitudes will be of interest. 

We see from (3) that the shape of e ± (x) in guide II is given by 

/(/3 ± , x) = [(pi/p 2 ) sin p 2 x + cos V&l\fi±- (23) 

Since 

/(/9 +I x) = f(0 2 , x) + (2A + 8) || 



s(fi^ X ) ± m,x)-i^ 



3i 



the shapes of both modes differ just slightly from the unperturbed 
shape /(/3 2 , x) ; furthermore, the shifts for the two modes are unequal 
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and are in opposite directions. The unperturbed shape can be deter- 
mined with the aid of (9). We find 

H? = tan \w 2 p 2 + l (tan- 1 ^ - tan" 1 & )] m tan U 

for even-numbered modes and 



P- 1 = - cot U 



for odd-numbered modes. Thus, 

f(0 2 , x) = tan U sin p 2 x -f- cos p 2 x 
= sec U cos (p 2 x — U) 
(P? + vlY 



Pi 



cos 



p 2 x — iv 2 p 2 



.Iftan-i^-tan-^)] (24) 
2\ Pi Pi)\ 

for even-numbered modes and, similarly, 

/(/ 3 2) X ) = 6* + pD* sin j^ _ Wap2 

-^tan-^-tan-^ 3 )] (25) 

for odd-numbered modes. The mode shapes in guide IV can be deter- 
mined in an analogous manner, with perturbations performed about 
/3 4 instead of jS 2 - We leave the details to the reader. 

The above results are not surprising. If the double waveguide 
system has a mode with propagation constant 0+ or /3_ which is close 
to the propagation constants /3 2 and 4 of modes that can travel in the 
individual isolated waveguides, then we would indeed expect the shape 
of that double waveguide mode to deviate only slightly in each wave- 
guide from the shape of the mode that could propagate in the isolated 
waveguide. 

The amplitudes of e ± (x) in guides II and IV prove to be more inter- 
esting. Let the arbitrary amplitude A in (3) be written as A ± for 
e±(x), and let B% and B*l denote the amplitudes of the modes in guides 
II and IV, respectively. Then (3) and (24) or (25) show that the mode 
amplitudes in guide II are given (to our order of approximation) by 

B u = A CgL±ifi!. (26) 

P2 
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In guide IV, the amplitudes are, by (3), 

B 1 ! = A ± \c 2 C 3 (l + £ T 2 ) (1 - XT 3 ) {vl + gS! 1 • (27) 

Care must be taken in the evaluation of C 3 (l — .XT3) at 0±. Since 
X(j8+) < 1, F(/3 + ) < 1, and by (8), 

T 3 = X+ J = tanh (tanh" 1 X + tanh" 1 Y), 

1 T A I 

we have 

1 + XF 
C 3 = cosh (tanh -1 X + tanh -1 Y) = (1 _ x 2 )Hl — P0> ' 

so that 

We find in a similar manner that 

C 8 (l-XT,)k= -[ (2A + j )r ,] • 
Thus the mode amplitudes in guide IV are 



(28) 



We observe that the mode amplitudes will have the same signs in one 
waveguide and the opposite signs in the other. Thus, e+(x) might be 
termed quasi-even and e-(x) quasi-odd. More startling, however, is 
the realization that, if A > 0, then the ratio \B^/BU\ = A+/A- may 
be quite different from the ratio IB^/B 1 !] = (A + /A_)[(2A + «)/«]. 
As we see in the next section, this will have serious implications when 
we consider the double waveguide system as a device for transferring 
power. 

For future reference, we write e±(x) for a system consisting of two 
symmetric, but not necessarily identical, waveguides {Ki = K 3 = K5). 
In this instance, (9) and (10) imply that 

C£L + (pi/pt)T x l = 1, 
C 4 [l + (pi/p*)TQ = 1, 
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so that we have by (3) and our previous analysis 

e+(x) = A + exp pix x < 

= A+ (p * + pi) * cos p 2 (x -w 2 ) < x < 2w 2 

= A+[— Xsinhpi(a; — 2w 2 ) + coshpi(x — 2w 2 )]|/3 + 

2u>2 < x < 2{w 2 + w s ) 

. r (2A + 8) X'V(pl + pl)i / o o \ 

= A+ - ^ — - — — *-=- cos p 4 (x — 2w 2 — 2w z — Wi) 

L 5 / J P4 

2(w 2 + w 3 ) < x < 2{w 2 + w 3 + w t ) 

= A+ \ (2A - + S) y, V exp Pl (2w 2 + 2w 3 + 2u> 4 - x) 

2(w 2 + w 3 + w^ < x, (29) 
e-(x) = A- exp p\x x < 

= ^_ ( P' ± P»)* cos p 2 (x -w 2 ) <x < 2w 2 

= A-{_ — X sinhpi(a: — 2w 2 ) + coshpx(x — 2iy 2 )]|/3_ 

2w 2 < x < 2(w> 2 + W3) 

2(w 2 + W3) < x < 2(w 2 + w 3 + W4) 

(2A + a)F J exp Pl( - 2w2 + 2 ™ 3 + 2m " x) 

2(w 2 + w 3 + wO < a:. (30) 

V. BEAT LENGTH AND POWER TRANSFER 

Suppose the two modes E ± = e±(x) exp i(ut — /3±z) travel down the 
double waveguide device. Since they travel at different phase velocities, 
the quasi-even and quasi-odd modes will alternately add and (par- 
tially) cancel in each waveguide. Hence, power is transferred between 
the two waveguides. 

The beat length L over which this transfer takes place is given by 

L " 2(5 + A) 25[1 + (A/5) 2 ]* ^ ' 

Note that, if the degenerate shift 5 is fixed, then as the mismatching A 
increases, the beat length L decreases. We can conceive of ways to 
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tune the double waveguide device and thus to change the beat length. 
This might be useful for light switching or modulation. 

It is important to learn just how much power can be transferred 
in the waveguide system. Suppose, for definiteness, that we excite just 
one waveguide at z = (say, guide IV), with the intent of transferring 
power to guide II via the mode-coupling mechanism. If guides II and 
IV have degenerate modes (A = 0), then as the modes travel down 
the waveguide system, they will alternately add and then cancel (to 
order 8 2 ) in each waveguide, with addition occurring in one wave- 
guide when at the same position cancellation occurs in the other. If 
guides II and IV have nondegenerate modes (A > 0), however, then 
complete cancellation cannot take place in both waveguides: by (27) 
and (28), we see that if the amplitudes of e+(x) and e_(x) are adjusted 
so that the modes cancel at z = in guide II, then they will never 
cancel fully in guide IV. 

From a practical point of view, a parameter which is likely to be of 
interest in this matter is the fraction of the total power introduced into 
the system which can be transferred into guide II. If the modes are 
poorly confined, an appreciable fraction of the power carried by a 
waveguide may actually be outside the high dielectric guiding region. 
If the reader is interested in the fraction of the power which can be 
transferred not only to the guiding region of guide II, but also to its 
vicinity, we would need a power transfer ratio G to bo defined by 



G = 



= j°_ [>+(*) + e.(x)Jdz J ' J' [_e\{x) + e 2 _(x)]dx, 

where a is some number between 2io 2 and 2(iy 2 + w 3 ) which defines 
the "boundary" between guides II and IV. The numerator of this 
expression, then, is proportional to the power carried by the entire 
guide II. 

Unfortunately, for a general asymmetric waveguide system, it is 
not at all clear how to define the position of the "boundary" between 
the two waveguides. If the system consists of two symmetric wave- 
guides which are nearly identical (except for a small deviation if the 
modes are slightly mismatched), then it seems clear that the boundary 
should be midway between the two dielectric regions, i.e., at 
a = 2w 2 + w 3 . By using (29) and (30), we find in this instance that 
we have to first order 

G = [1 + (A/8)*3-K 

Thus for perfectly matched waveguides (A = 0), the power transfer 
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is complete, to first order. As the mismatching increases, the power 
transfer ratio decreases rapidly. 

Complete power transfer (to first order) is a direct consequence of 
assumption (17), which implies little overlap between the field as- 
sociated with guide II and that associated with guide IV. A higher 
order perturbation analysis would show that in fact there is some field 
overlap and that, even if A = 0, the power transfer is not complete. 
As the waveguide separation increases, there would be less field overlap 
and the power transfer would be more nearly complete. 

For a general asymmetric waveguide system, we might define the 
"boundary" between the two waveguides to be, say, at the position 
where the "quasi-even" field attains its minimum. Such a definition 
can be cumbersome to apply mathematically. In general, though, we 
would expect results similar to those obtained for the symmetric 
system. If the modes are degenerate and one waveguide is excited, then 
virtually all the power can be transferred to the vicinity of the other 
waveguide. The power transfer ratio decreases as the mismatching 
increases. 

It will be instructive to introduce a second power transfer ratio F, 
which can be defined precisely. It will be the fraction of the total power 
introduced into the system which can be transferred into the high 
dielectric region of guide II, the waveguide which was originally un- 
excited. If terms of order A//3 are neglected, this power transfer ratio 
is defined by 

F = J [>+(*) + e-{x)Jdx/ j_ [e 2 +(x) + e 2 -{x)~]dx, (32) 

where the mode amplitudes A+ and A~ are equal. 

If the modes are poorly confined in guide II, the power transfer 
ratio F may be considerably less than unity even if the waveguides 
are perfectly matched (A = 0). The definition of F is concerned only 
with the power which can be transferred into the high dielectric region 
of guide II; hence, F depends upon the confinement factor of the 
waveguide (to be defined below) as well as upon the amount of mis- 
matching A. 

Evaluation of (32) can be very messy for the general case of two 
asymmetric waveguides. We simplify the subsequent analysis and yet 
retain its essential flavor by assuming that the double waveguide 
structure is composed of two symmetric, but not necessarily identical, 
waveguides (ifi = K 3 = Kb). The modes e ± (x) are then given by (29) 
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and (30). For the numerator of F, we find 

/ [e + (.c) + e-(x)Jdx = 4A 2 + KPl ^ P2) / cos 2 p 2 (.r - iv 2 )dx 
Jo Pi J 

= 4A 2 + (p * + 3 Pi) [p 2 iu 2 + * sin 2p 2 u; 2 ]. 
P2 

But since (9) implies that 

2pip 2 



sin 2/} 2 i<; 2 = 



v\ + vl 

for a symmetric waveguide, we have 

Similar procedures can be used in the evaluation of the denominator 
of F. Special care must be taken in evaluating the integral in the 
interval 2zt> 2 < x < 2(w 2 + io 3 ), where e+(x) and e_(.c) have the same 
functional description, but the function is evaluated at /3+ and /3_, 
respectively. In the resulting analysis, C 3 and £3 must be evaluated at 
/3 ± . Procedures similar to those used following (27) are helpful. The 
power transfer ratio turns out to be 

p = 4(P„/p 2 ) 

" 2[(Pn/p,) + (1/pi)] + 2[(P, V /P4) + (1/Pi)](l + 2AV8 2 )X'/r ' 

where 

p.™(^I) + £, 

P„_^(rf + P?) + *. (33) 

The expression for F can be simplified. By using (19), (20), and 
(33), we find that 



p 2 pi F' \ P* Pi J 



so that 



' [(Pn/p.) + d/Pi)][l + (A 2 /5 2 )] ^ ; 

Both (34) and (35) have interesting physical interpretations. In 
order to discuss them, we make a brief digression. Suppose that, 
instead of the double waveguide system, we just have an isolated 
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guide II, which for generality we assume is not necessarily symmetric. 
(In the notation of Section II, we have K 3 = K t = K 5 .) If the electric 
field is given by e(x) exp i(wt - /3 2 z), then the power carried by the 
entire waveguide is proportional to 



= / e 2 {x)do, 

J—oo 



and the fraction of the total power which is confined to the high 
dielectric region is given by 

C = e 2 (x)dx I / e 2 (x)dx. 

It is not difficult to show that 

p = P?±_pi [2 + x + ( Pl / Pi )] 

where e(x) was assumed to have unit amplitude at x = 0. Upon 
comparing this expression with (19), we find that 

Y , _ 2^o (gj + pi) p 
Ps (p? + Pi) ' 

so that X' is related in a simple manner to the power carried by the 
waveguide with which it is associated. Now if the isolated guide II 
happens to be symmetric, then it is also true that 

Vi Pi 

X' = -^°P. 
Pi 

Since a similar relationship holds for an isolated guide IV, (34) follows. 
The physical implication of (34) is that, if a single waveguide is ex- 
cited in a matched double waveguide system (A = 0), then the power 
is distributed between the two modes in such a manner that 



f M e\(x)dx = f e 2 -(x)dx. 

J —at} J —oo 



If an isolated guide II is symmetric, then the confinement factor C 
can be shown to be given by 

n - jWPj (9A\ 

L ~ (Pn/p.) + d/Pi) V 
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Hence, the power transfer ratio F is simply given by 

F = C[l+ (A/5) 2 ]- 1 . (37) 

Note that F is completely independent of the parameters of guide IV, 
as well as of w 3 because of assumption (17). 

We remark once more that these results are first-order approxima- 
tions. If a higher order perturbation analysis were undertaken, it 
would show that F is also dependent upon the amount of field overlap 
between the two waveguides. 

It is important to note that the amount of mismatching A can have 
a significant effect on the power transfer ratio : as A increases, F de- 
creases rapidly. We might remark that, while the beat length L 
depends upon 5 and the ratio A/5, the power transfer ratio F depends 
only upon A/5. Thus, by proper device design, it may be possible to 
adjust 5 and A/5 to get both appreciable power transfer and a desirable 
beat length. 

VI. A NUMERICAL EXAMPLE 

To illustrate our results, let us give an example, using parameters 
which could be realized in a GaAs — AUGai-zAs heterostructure. We 
consider two waveguides which are each symmetric, but which have 
different widths and dielectric step heights. Our intent is to excite one 
waveguide and then to transfer power into the other by means of mode 
coupling. The parameters of one waveguide, say guide II, are taken to 
be fixed. The width of the other waveguide is considered a variable. 
For any given width 2w\, we adjust the dielectric height K\ so that 
the propagation constants /3o and fit for the zeroth order TE modes of 
the two waveguides match. 

We learn how the degenerate shift 5 varies with the spacing 2w; 3 
between the two waveguides and with the width 2io i of guide IV. We 
look at the beat length L and the power transfer ratio F for an idea of 
how much mismatching of the propagation constants /3 2 and 04 can be 
tolerated. We then discuss the amount of mismatching 2A = | /3 4 — 02 1 
which might occur in a practical situation and show how tuning can 
reduce this. 

To be specific, suppose that 

K 2 = 11.868, 
Ki = 0.9K 2 = 10.681, 
w 2 = 0.1 /im, 
A; = 5.4494 X 10 4 cm" 1 . 
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The dielectric constant K 2 corresponds to an index of refraction 
K\ = 3.445 of GaAs at a wavelength X = 2tt k~ l = 1.153 Mm. The 
halfwidth w 4 of guide IV will be assumed to vary between 0.1 ^m and 

0.8 f*m. 

For the zeroth order TE mode of a symmetric waveguide, we find 

that (9) reduces to 

tan p 2 w 2 = P1/P2, 

Pl = 031 - WKi)*, 

p 2 = (k>K 2 - /3f)». (38) 

The above parameters for guide II give us 

fr = 1.8049 X 10 6 cm" 1 , 
p x = 2.9311 X 10 4 cm" 1 , 
p 2 = 5.1631 X 10 4 cm- 1 . 

Equation (38) holds for guide IV if each subscript 2 is replaced by a 4. 
For any given value of w i} we adjust K4 so that /3 4 is the same as /3 2 . 
Some values of K A and p 4 are given in Table I. 

The parameters for K h K 2 , w 2 and any given pair K it u> 4 then define 
two waveguides which have a degenerate mode at the given wave- 
length X. The degenerate shift 5 is given by (21) for symmetric wave- 
guides : 

Pip 2 p 4 (l - tanh WzPi) 



5 = 



20 2 [(p? + PIKP? + Pl)(l + P1W2MI + PiW4)J J 



Figure 2 shows 8 as a function of w 3 for various values of w 4 . We see 
that the coupling decreases rapidly as the waveguide separation in- 
creases or as the width of guide IV increases. The beat length L, then, 
will increase rapidly with w 3 or w t if A is small enough. If, however, the 

Table I — Values of K 4 (and hence p 4 ) needed to match /? 4 for the 

zeroth order TE mode in guide IV (of halfwidth w 4 ) 

to the /? 2 for the corresponding mode in guide II 



w t (mitx) 


K t 


Pi (10* cm" 1 ) 


0.1 


11.868 


5.1631 




0.2 


11.381 


3.4917 




0.3 


11.222 


2.7340 




0.4 


11.145 


2.2763 




0.5 


11.100 


1.9619 




0.6 


11.071 


1.7295 




0.7 


11.051 


1.5494 




0.8 


11.037 


1.4048 
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w 3 IN Jim 

Fig. 2 — The degenerate shift 8 as a function of half the distance between the two 
waveguides w 3 , with the halfwidth w t of guide IV as a parameter. 



waveguides cannot be fabricated to match as closely as desired, the 
story can be different. If we take, quite arbitrarily, A = 100 cm -1 , we 
would find for example from (31) and the data in Fig. 2 that if 
Wi = 0.6 nm, 



L = 



0.113^ 
0.137 
0.150 
0.157 



> mm 



if 



103 = i 



0.41 
0.5 

0.6 



Mm (A = 100 cm" 1 ) 



while if A = 0, 



f 0.162 1 
L - i 0.279 I mm if 103 = 
0.490 




ixm 



[A = 0). 
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Thus, L is reduced significantly and changes less rapidly with w 3 if A 
is large enough. Similar results are found if w a is fixed and w 4 varies. 

Next, suppose we wish to excite one waveguide and to transfer 
power to the other one. If the propagation constants for the two 
waveguides were perfectly matched (A = 0), then by (35) to (37) 
the power transfer ratio F would be the confinement factor C which 
is plotted in Fig. 3. The upper curve is used if guide II is excited and 
power is transferred to guide IV; the lower curve is used if guide IV is 
originally excited and power is transferred to guide II. The power 
transfer ratio is larger if power is transferred from a narrow guide II 
to a wide guide IV than if the power is being transferred the other 
way, since the power is more tightly confined within the guiding region 
in the wider waveguide. 

Lest the reader become confused, we recall that F is defined as the 
fraction of the total power introduced into the system which can be 
transferred into the high dielectric region of the guide which was 
originally unexcited. F does not concern itself with how much power 
in the high dielectric region of one guide can be transferred into the 
high dielectric region of the other guide. 




0.4 0.6 

w 4 IN pm 



Fig. 3 — The confinement factor C for degenerate modes. Guide II is fixed, and the 
half width w t of guide IV varies. 
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Again, the mismatching A can have a significant effect. If, for ex- 
ample, we assume that w 4 = 0.6 /im and transfer power from guide II 
to guide IV, we find from (35) (evaluated with the parameters for 
guide IV) that if A = 100 cm" 1 , 



F = 



0.66] 




ro.4 


0.44 L 


if Wa = ■ 


0.5 


0.22 




0.6 



Mm. 



(The maximum value of F is 0.91 for A = 0.) 

If we know that the fabrication procedure will likely make A of 
significant size, then the only way (for given waveguide parameters) to 
get good power transfer is to make 8 large enough. A trade-off thus 
must be made between good power transfer and a long beat length. 

Fortunately, tuning can be a viable alternative to making such a 
trade-off for badly matched waveguides. To tune a device which has 
already been fabricated, we would need to alter one or more slab 
widths or dielectric heights. 

Some possible methods of tuning are to change K 2 or K\ by altering 
the free carrier density or using the electro-optic effect, or to change the 
outer slab levels Ki or K 5 by diffusion or ion implantation. This can 
be achieved, in principle, by growing at least one waveguide with a 
small gradient in the slab width. Phase matching then can be achieved 
by lateral positioning of the light beam, which travels approximately 
perpendicular to the gradient of the slab width. 

We compute a possible value of A for a specific example and then see 
how much tuning is needed to reduce A to zero. Suppose that (for any 
given w t ) the double waveguide device is fabricated according to the 
specifications for matched propagation constants, but that there are 
slight errors in w 2 , Wi, K 2 — K h and i£ 4 — K\. Assume the following 
errors, which are probably reasonable if the device is fabricated by 
molecular beam epitaxy: 5 the ratio w 2 /wa is nearly constant, and w 2 
varies by ±0.02 nm; the ratio (K 2 — Ki)/(Ka — Ki) is nearly con- 
stant, and K 2 — Ki varies by ±0.10. We shall take K 2 to be fixed. 
Then the extreme cases would be given by w 2 = 0.12 (0.08) /zm, 
K\ = 10.781 (10.581), and 

w\ = (w i /lV 2 )W2, 

Ki - Kl ~\K 2 - Kj Kl + \KT=- Kj Ki ' 

where the primed parameters refer to the values in the fabricated 
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Fig. 4 — The mismatching A as a function of the ratio of the guide widths. 

device. We just treat the extreme case for which w 2 = 0.12 /xm and 
K\ = 10.781, since the other gives changes of comparable size. 

We find from (38) that, with these fabrication errors, /3 2 = 1.8156 
X 10 5 cm -1 and /3 4 varies from 1.8156 X 10 6 cm -1 for w 4 = w' 2 to 
1.8119 X10 6 cm -1 for w 4 = 8w 2 - The resulting values of A = §|/3 4 -/3 2 | are 
plotted in Fig. 4. Although both fa and 4 have changed significantly 
from the value 1.8049 X 10 5 cm -1 for which the device was designed, 
they both change in the same direction, so A reflects a less radical 
change. 

In tuning the system, suppose we consider guide IV, and hence 4 , 
as fixed for any given wi ; we shall alter either the dielectric height or 
the symmetry of guide II to adjust /3 2 . 

If we lower K 2 to make fit match 4| we find by (38) that the altered 
values K 2 are those given in Table II. The change in K 2 is thus less 
than 0.4 percent. Such a change is feasible with free carrier injection 

Table II — Values of K 2 needed to tune the mismatched 
waveguide system 



w'Jw'% 


K\ 


% change 


1 


11.868 





2 


11.846 


0.19 


3 


11.840 


0.24 


4 


11.831 


0.31 


5 


11.829 


0.33 


6 


11.825 


0.36 


7 


11.824 


0.37 


8 


11.823 


0.38 
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Table III — Values of Ki needed to tune the mismatched 
waveguide system 



w\/w' 2 


K[ 


% change 


1 


10.781 





2 


10.716 


0.6 


3 


10.697 


0.8 


4 


10.670 


1.0 


5 


10.663 


1.1 


6 


10.657 


1.2 


7 


10.647 


1.2 


8 


10.643 


1.3 



and is about an order of magnitude larger than can be handled com- 
pletely by the electro-optic effect. 

Another possible method of tuning would be to make guide II 
asymmetric. If we keep K 2 = 11.868, K3 = 10.781, w 2 = 0.12 /xm, and 
alter K\ to match /3 2 to /3,», we use (9) in the form 

■pi tan 2W2P2 — Pi 
1 -I- (P3/P2) tan 2w' 2 -p 2 

to get pi and then use (4) to find K\. The altered values K\ are given in 
Table III. The change is thus no more than 1.3 percent. This can be 
handled by diffusion or ion implantation. 
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